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AhSTKACT 


The atruoturo of cometary nuclei Ir a key Issue tn solar system 
cosmogony. The lntern.il structure of a comet could he determined with 
a spacecraft borne electromagnetic sounder. A dielectric profile of 
the comet could be produced In direct analogy with terrestrial glacier 
<ind Ice sheet sounding experiments. This profile would allow the detec- 
tion of a rocky core or lea layers if they exist, Just as layers In the lee 
and the bedrock .‘nterface have been clearly observed through the Greenland 
Ice sheet. It would also provide a gross estimate of the amount of dust 
In the Icy region. Models for the response of the nucleus and comet.«i y plasma 
to electromagnetic sounding are developed and used to derive experimental 
parameters. A point system design was completed. Preliminary engineering 
study t'esults Indicate that the sounder Is well within the bounds of current 
sp.ice technology. 
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SECTION 1 
INTRODUCTION 


ORIGINAL PAGE IS 
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The origin and structure of comets is a key Issue In solar system 
cosmogony. One theory suggests that comets are In fact nearly unaltered 
samples of the primordial material from which the sun and planets accreted, 
and hence, beer information about the conditions extant <4.5 billion 
years ago when the solar system was formed. Alternatively, it is possible 
that comets contain material from the interstellar medium and may predate 
the solar system ( Report of the Comet Halley Scienc e ulorking Group. 

1977). Comets are known to contain organic molecules (Ip and Mendis, 

1976); these molecules may have played important roles in forming the 
very early environments of all the planets. Additionally, what relation 
comets have to asteroids and 'eteorites, aside from well-known orbital 
similarities, is unknown. It is possible that asteroids and meteor 
streams are defunct comets, long since stripped of whatever volatile 
ingredients give comets their spectacular appearance. Since current 
tl eories of the solar system's origin are based in part on data from 
these asteroids, what relation they have to comets is a question of 
paramount importance. 

Current knowledge of comets is quite scarce, and is based almost 
solely on spectroscopic analysis of light emanating from the coma and 
tail. These data have allowed tentative identification ol various 
organic-free radicals (Ip and Kendls, 1976; Whipple and Heubner, 1976) 
and considerable speculation about the molecular origins thereof. 

Some theoretical models covering physical ionization processes in the 
coma and the interaction of the coma with the so’ar wind have be-'n 
attempted, but there is no direct information on the size, structure 
and composition of comets' nuclei and the cometary plasmas, and hence 
these models are somewhat speculative. 

The C'^metary model with the widest endorsement is undoubtedly 
that of Whlp-'le (Whipple, 1950). A comet is postulated to con.-^ist 
of a large nucleus made of rocks with substantial amounts of accreted 
ice, typically several kilometers across. Tnere may be lar«o ouantlties 
of impurities, dust and chemical matter imbedded in the ice. In thl- 
model , when a comet approaches the sun, the nuclear ice starts to sublime 
and ionize, thereby dispersing dust and molecules, which are blown 
by the solar wind to form the cometary tails. 

Direct missions to comets are needed to finally solve some of 
the key problems in cometary structure rese'.rch (EepLOTt .Qf Jhe. .Conict 
halley :jcience Workin^p Gr ou p . 1977) such as: Verify that a solid nucleus 

exists at all, and if it does exist, what is its size? What is the 
nature of the nuclear ice? Does the nucleus have on# large rocky cor« , 
a collection of smaller rocks .»id du.st embedded in ice, or some comt)ina- 
tion of the two? If there is a core, what is its size and share? 

The objective of this report is to discuss the potential of a 
spaceborne electromagnetic sounltr to answer some of these questions 
and to evaluate the feasibility oi such a sensor . 
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fRCTION II 


ELECTHOMAGNETIC SOUNDING OF THE NUCLEUS 


Deep electromagnetic sounding of the nucleus from a nearby location 
(i.e., less than a few tens of kilometers) could unambiguously resolve some 
of the key questions about nuclear structure; in parti sular, the presence 
or absence of a rocky core and its size, the presence of layering in the 
ice, and the spatial distribution and amount of rockj and dust in the icy 
region of the nucleus. 

Electromagnetic sounding of a solid nucleus is, at least in 
principle, very simple. An antenna launches a short-pulsed electro- 
magnetic wave aimed at the nucleus. The wave propagates through space 
and strikes the nuclear surface. Part of the wave is transmitted, and 
part is reflected. The transmitted wave continues to propagate through 
the nucleus and is reflected by any macroscopic change in the dielectric 
constant. Thus, a rock-ice interface or a large center of Impurities 
will cause some reflection. Meanwhile, the antenna is switched into a 
receiving mode and will detect any radiation that is reflected back. 

The first return will come from the surface Interface; later returns 
will come from scattering centers and interfaces inside the nucleus, 
and there may be some reflected radiation from the rear surface interface 
(Figure 2-1). By probing from different locations around the nucleus, 
a three-dimensional "map" of the internal structure may be deduced. 

Sounding techniques have been developed for terrestrial glaciers 
and ice sheet mapping (Gudmandsen, 1975, 1976; Robin et al . , 1969; 

Hockstra and Capillino, 1971; Evans, 1965). Sounding of Greenland and 
Antarctic ice sheets, down to the bedrock through many kilometers of ice, 
has been achieved with airborne electromagnetic sounders (Figure 2-2) 
using a reasonably low amount of power which is compatible with spacecraft 
capabilities. Numerous layers have been observed throughout these ice 
sheets. Water ice is nearly transparent for electromagnetic radiation 
at frequencies between 0.1 and 100 MHz (Figure 2-3). Waves launched 
at these frequencies will propagate with little loss through the ice 
and reflect from rocks, bubbles, dislocations and finally from the 
underlying bedrock. The time separation between the reflections from 
the surface and bedrock interfaces allows the determination of tne 
ice thickness assuming that the real part of the ice dielectric constant 
is re?iSonably well known. As the sensor platform moves above the surface 
of ti.-j ice, a complete profile of the internal structure is obtained. 

The extension of the above method to the sounding of a comet nucleus 
is straightforward; the principal component of nuclear material is expected 
to be water ice (Ip and Mendis, 1976), so the above conclusions about pro- 
pagation through glacial ice sheets hold reasonably well. The effect of dust 
in the cometary ice will be discussed later in some detail. The radio sounder 
could measure the size of the nucleus and determine whether or not it has a 
core. Sounding could determine the approximate value of the dielectric con- 
stants of the surficial and internal material. By probing the nucleus from 
many different directions, the location of small rocks, dust patches and other 
scattering centers might be determined in direct analogy with glacier sounding. 
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Figure 2-1. Specular Return from the Nucleus. Electromagnetic 
waves launched from the antenna are reflected and 
transmitted at the surface of the nucleus (front 
interface). The transmitted energy propagates through 
the nucleus and is reflected at its back surface. 
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Figure 2-2. Internal Structure of Greenland Ice Sheet. The flight 

track of a plane carrying an electromagnetic sounder corresponds 
to the horizontal axle. Depth is on the vertical axis. 

Layers and the bedrock interface are clearly seen. Sounding 
to a depth of more than 3 km has been reported. The sounder 
frequency is 60 MHz. (Courtesy of the Electromagnetic 
Institute, Technical University of Denmark.) 
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2-3. Attenuation of Electromagnetic Wave? Ir Ice. 

Absorption per 100 n if, giv< n as a functjon < * 
Irequency, with temperature as a parameter (li * 
Smith and Evans, 1906). 
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SECTION III 


^REQ0ENCy SELECTION 


Th* aaplitud* of th« signal backscatterad froa the nucleus will 
alaost inevitably be reduced aany orders of aagnitude froa the aaplitude 
of the transaitted radiation. Coapounded with liaited power typical 
of spacecraft sensors, this aeans the sounder aust be designed to maxiaize 
the returned signal. It is particularly iaportant that the antenna 
has a high directivity gain and that the transaitted radiation is nowhere 
evanescent in the coaetary plasaa. 

This latter requlreaent iaaediately precludes sounding the nucleus 
with low frequency radiation. In particular, electroaagnetic waves 
will not propagate through the plasaa at frequencies below the plasma 
frequency, upt 



(3-1) 


where 


n a electron density 
m a electron mass 


e a 'Electronic charge 

In a cometary plasma, the maximum electron density will occur 
at perihelion, when the nuclear ice is sublimated at ». maximum rate. 
According to Ip and Mendis (1976), a typical peak electron density 
at a 1 All perihelion for a 5-km nucleus on the sunward side is 10^ cm~3. 
Since observed plasma isophotes are approximately elliptical (Feldman, 
1977), it is reasonable to conclude that the peak electron density 
is not any lower on the dark side of the nucleus. 

The plasma frequency corresponding to an elect"on density of 
10^ cm"3 Is approximately 10*^ Hz (i.e., 10 MHz). This may be regarded 
then as an overall lower limit on the sounding frequency from any sounding 
position around the nucleus. 


Higher sounding frequencies are advantageous because they require 
smaller antennas; however, care must be taken that the attenuation 
of the wave in the ice does not become too severe (Figure 2>3)i and 
that scattering by rocks with a size of a fraction of a meter or less 
does not attenuate the wave excessively. The above factors indicate 
that a sounding frequency range of 30 to 100 MHz is the most desirable. 


The impurities in the nuclear ice are most probably ammonia, 
methane, carbon dioxide, and dust (Feldman, 1977; Whipple and Heubner, 
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Wt) . It is also possible that laore axotlo molecules are present, 
but •>nly In trace amounts, too small to affect the propatcatlnc electro- 
' a>'nrtlo wave. Ammonia, methane, and carbon dioxide, even if they 
irr quite abundant, will only marginally change the refractive index 
cl t^ater ice (Mestphal, 1977) a.id may be Ignored in preliminary cal- 
culations. It is a reasonable assumption that no compounds with anomalous 
dielectric behavior are formed from these impurities. Dust, however, 

.nay increase the dielectric constant in the nucleus substantially; 
there may be as much as a factor of two difference between the real 
•dielectric constant of ice and that of dirty ice with silicates or 
"rocks." The effect of dust concentration on the electromagnetic waves 
is mo<1<'led in the next section. 
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SECTION IV 

HOOLLS PHCOICTING THE RETUKNED SIGNAL 


A number of phenomena would reault in a returned echo ac the 
incident electromagnetic wave propagates through the nucleus: 

(1) The wave is reflected at the front surface of the nucleus, 

at any interface (such as layers in the ice or on the surface 
of a rocky core) throughout the nucleus, and at the back 
surface of the nucleus (see Figures 2-1 and k-la). As it 
propagates through the nucleus, the wave is attenuated 
by absorption in the ice, the dust and the impurities 
and scattered from reeky fragments. 

(2) The wave can be refracted and returned back to the antenna. 
This is jomewhat similar to the "halo" effect and Luneberg 
lens effect (see Figure ^4-1b). 

(3) The wave can be backscattered from any roughness on the 
surface of the nucleus (see Figure 4-1c) . 

Other phenomena can also occur, such as surface waves which would 
circumnavigate the nucleus and re-radiate back toward the antenna 
(Figure 4-1d). However, in this preliminary report wc will only consider, 
in a simplified fashion, the first three phenomena listed above, which 
seem to be the most important. 


A. HEFLECTION A.'.L.' AbSURPTION 

Backward reflection occurs when a wave is normally incident on 
an interface where there is a change in the dielectric constant. The 
exact calculation of the backward reflection from a layered dielectric 
sphere Is somewhat involved. However, a planar layered medium can 
be used (see Figure 4-2) to obtain a first-order estimate of the echos 
from the different interfaces in the nucleus. 

In the case of a homogeneous nucleus (Figure 4-2a) , the power 
Pp of the echo reflected by the front surface is equal to 


PtGAri 


r 5 

I6n(d - a)^ 


(4-1) 


where 

P^ s transmitted power 
G : antenna gain 
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Figure 4-1. Phenomena Causing an Echo, (a) Reflection and 

scattering from Internal inhomogeneities including 
layers, large core or relative! i small-scale 
inhomogeneities, (b) Uniquely refracted modes where 
the wave path is refracted inside the nucleus and 
returned straight back to the source, (c) Back- 
scattering from the surface roughness, (d) Surface 
waves which would propagate all around the nucleus. 
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Figure 4-2. Nuclear Models and Corresponding Planar Models. 

(a) Homogeneous nucleus of ice (including dirty 
ice), (b) Nucleus with a rocky core, (c) Layered 
nucleus, (d) Nucleus with a rocky crust and an 
icy core. 
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A X intenna effective receiving area 
d X distance to the center of the nucleus 
a s radius of the nucleus 

r^ I Fresnel reflection coefficient of the surface 





The power of the echo fron the ‘?a'k surface is equal to 


I 



GAr2Vl2 



where 

I X transmittivity of the first interface 
V = roundtrip attenuation factor for the nucleus material 
r 2 = reflectivity of the second interface 
Uf particular interest is the power ratio of the two echos: 


(4-2) 


I 



d - a 


.1 ♦ 


if-') 


Ta 


(4-3) 


In the above expression, we took = r 2 - This ratio is plotted in 
Figure 4-3 for different values of d/a and for * = 3.2 which is the 
dielectric constant of water ice. This ratio represents the strength 
of the echo fron the back sHe of tne nucleus versus tne echo fron 
the front siie. In other t.*rras, it is the capability of seeing all 
through a homogeneous nucleus. To Illustrate, let us assume that d/a = 
4, a = 2.5 km and the absorption in the loe is 3 dB/kra. Then we have 


^r 

— = -33 (IB 
f'r 
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Hgure 4-3. Echoes from Eront and Rear Nuclear Surfaces. Ratio of power 

return from the back interface (Pp) and front interface of the 
nucleus (Pp) versus the ratio of the distance of the spacecraft 
from nucleus center (d) to nuclear radius 'a). Hie roundtrip 
absorption in the nucleus is a. 


Thus, if we impose that the backside echo miwt have a signal -to-noi.ie ratio 
(Shh) of at least 13 dP, then the front side echo must have an SNR of 46 dP. 

In the case of a rocky core (Hgure 4-2b) of radius Rg, the power 
ratio of the echoes from the nucleus surface and the rocky core surface 
is simply equal to 



«here r^ is the reflectivity at the icc-rock interface. 

This ratio is plotted in Figure 4-4 for a rocky core of dielectric 
constant equal to 7 which is typical of silicate materials. The echoes 
IVom the back surface of the core and nucleus are expected to be very 
weak because of the absorption which occurs as the wave propagates through 
the rocky core itself. 

If we assume that the core has a 1-km radius, then using the 
same parameters as above, we find 





-12 dB 
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Figure U-4. First Echoes from Nuclear Surface and Core. 

Itatio of power return from the core-ice Interface 
(Pp) and the nucleus surface (Pp). The roundtrip 
absorption in the ice crust is u . 


Thus, if we inpose a minimum SNR of 13 db for the echo from the core, 
then the front side echo must have an SNR of 25 db. 

The .SNR imposed on the front echo defines the parameters and 
performance of the sensor. In Table U-1 we give the parameters of 
a point design sensor which is well within the present state-of-the-art 
of space technology. To illustrate the performance of such a sensor, 
let us assume d s 10 km. Then 


SNR s = 80 db 

kTB 


where Pp is given in Equation (4-1). This gives a very high SNR, which means 
that sounding could be easily accomplished from even farther distances. 

This will be discussed later. 

It is possible that each time a comet passes through perihelion, 
only the volatiles in the icy mantle are stripped off, and that dust 
and rocks are left behind. Particularly in an older comet, this could lead 
to a "rocky shell" effect (Figure 4-2d), where a high concentration of rocks 
and dust remain near the surface of the nucleus with an icy core at the center. 
Using the sensor parameters in Table 4-1, the SNR predicted for the echoes 
from an equivalent slab model (Figure 4-2d) is shown in Figure 4-5. The 
SNR was calculated for the case where the spacecraft is 10 km away from the 
surface of t^e nucleus. It is clear that even with a 1-km thick crust, 
complete sounding can still be achieved from a 10-km range. 
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Table <(-1. System Parameters 


Parameter 

Value 

Frequency 

f 

S 

*40 MHz 

Wavelength 

A 

S 

7.5 m 

Peak Power Transmitted 

Pt 

B 

500 W 

Pulse Length 

5 MS 


bandwidth 

B 

S 

2 MHz 

Range Resolution 

75 

m 

in free space 


**5 

ffl 

in ice 

Ailse Repetition Frequency (PRF) 

10 



Average Transmitted Power 

0. 

1 Vi 

1 

Receiver Noise 

T 

s 

600 K 

Antenna 

Yagi 

(3 dipole elements, 


each 

1.9 m long); main 


lobe 

has azimuthal symmetry 

Antenna Gain 

G 

= 

6 db 


The most general case consists of a large number of layers with a 
continuously increasing dielectric constant toward the center of the 
nucleus (Figure 4-2c). This case can be easily simulated if the curvature 
effects of the incident wave and the nucleus are neglected. Figure U-2c 
illustrates one case wliere the dielectric constant varies linearly in 20 
layers from « : at the nucleus surface to c = ^4 at the nucleus center. 

Figure *4-6 shows the SNR for reflection from this layered slab racdel, whore 
the envelope of the peaks from the successive interfaces is plotted. 


b. fciFECTS Oh' DUST AND IMPURITIES 

To study the effects of impurities we calculated the change in 
the dielectric constant of ice when impurity with a different dielectric 
constant is added. The resulting composite dielectric constant is 
given by (Von Hlppel , 195*4) 
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Klgure 4-5. First and Second Echoes from the Nuclear Surface and 
Core. Slgnal-to-nolse ratio of the echoes from the 
lour Interfaces which correspond to a nucleus with 
an icy core and a rocky shell. The sounder is 10 km 
away from the nucleus. corresponds to the echo from 
the front surface. Ej is for the echo from the front 
rock-ice interface. E 3 is for the echo from the back 
ice-rock interface. Ejj is for the back surface. R(, 
is the core radius. 

log » s p log Ic ♦ (1 _ p) log (, (4-5) 

where 

t = composite dielectric constant 
p = impurity concentration by volume 
' imp = impurity dielectric constant 
’ ice - i®® dielectric constant 

Ttie composite dielectric constant for various Impurity dielectric 
constants is shown as a function of impurity concentration p in Figure 4-7. 

Audit iona 1 1 y , impurities affect the absorption properties of 
lee. The attenuation of the power of the wive in a lossy medium due 
‘o hayieigh absorption by tiny distributed .scattering centers is given 
t y (hattan, 1972 ) 


4-B 



‘•HIGINAL pagk is 
»»K IH»OR quality 



Ki(?ure 4-b. 31^nal-to-Noise Ratio for a 20-Layer Nucleus. The curve Is 

the envelope of the echoes from each Interface. The horirontal 
axis Is the time between the transmitted sii^nal and received 
echo. The sharp rises at 3** mo ind 08 ps correspond to the 
echoes from the front and bacK surfaces respectively. The soiimler 
is 10 km from the nucleus. 



Figure 4-7. Dielectric Constant of Ice with Impurities. The 

dielectric constant of ice is plotted as a funct ion 
of Impurity concentration (by volume) for different 
types of impurities. Curve A corresponds to = 4, 
tan fti = 10"3; curve B corresponds to = 5, tan <S 
10“2; curve C corresponds to = 8, tan b = 10”^; 
curve D corresponds to = 10, tan h = 10“\ 
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where 


Fp X power at distance r 
PpQ X initial power 

N X number of lossy scatterers/m^ 
r X distance into medium 


Og X 4nka^ Im 



1 



k X 2v/\, X X wavelength 
a X scatter size 

m X index of refraction (square root of dielectric constant) 
of scatterer 


The standard form for attenuation by a slightly lossy medium is 


exp 


|?ki ^1 


♦ i tan 6 


) * 


exp 


( i tan 6 \ 
'* — ) 


where tan 6 is the loss tangent of the medium. Equating the real parts 
of the exponents, we obtain the equivalent loss tangent 


tan 6 = -SnNa"^ Ira 


m^ - 1 
m^ ♦ 2 


We also have 


Ira 


m^- 1 
ra^ ♦ 2 


3c" 


3c' 


(c* 2)‘ 


(c * ♦ 2) 


2 ^imp 


where tan s < "/ i • 

Thus 


l8i • 

tan 6 = p tan 

(c* ♦ 2)" 


(4-7) 
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where p la the percentage of lapurltlea by volume. This expression 
gives the loss tangent due to Impurities. The total equivalent loss 
tangent also Includes the absorption In the Ice matrix. This gives 


Un H s 


I8i • 

( 1 - p) tan ^ice ♦ P ■ ■ * * * x ^imp 

(.' ♦ 2 )^ 


(*»-8) 


A plot of tan for various values of Impurity concentration Is given 
In Figure <t-8. For comparison, we give In Table 4-2 the values of tan h 
for a number of minerals encountered on Earth. 

To evaluate the effect of dust concentration on the capability 
of sounding the nucleus, let us consider the case of a homogeneous 
dirty Ice nucleus. We plotted In Figure 4-9 the SNR of the echo from 
the back of the nucleus (i.e., opposite side from the spacecraft) as 
a function of the loss tangent. If we assume a minimum SNR of 13 db, 
then the equivalent loss tangent must be lower than 1.1 x 10*3. This 
sets an upper limit on the amount of impurity which would still allow 
sounding of the entire nucleus (from Figure 4-8). 

If there are well-defined layers, the returned signal would have 
a number of well-defined echoes similar to the case of Earth ice sheet 
sounding (Figure 2-2). As the spacecraft moves around the nucleus, 
a three-dimensional model would then be obtained about the spatial 
extent of *.hese layers. 

If the nature of the nuclear material changes gradually toward 
the center, the returned signal would consist of a diffuse echo (as 
opposed to distinct (ulses). The shape of the echo could then be used 
to deduce information about the internal structure of the nucleus. 

This latter case would be somewhat harder to interpret exactly. 

In actuality, the analysis of the sounding data, which would 
be obtained during a comet rendezvous mission, is an inverse problem. 

The internal structure is unknown while the echo characteristics are 
measurable. The only other relevant information available would be 
the size of the nucleus (from photograpt.y and ranging), some idea about 
the composition at the surface from analysis of the ejected material 
(bid possibly the nuclear mass. 

If the nucleus is a homogeneous ice ball (Figure 4-2a) , the returned 
signal would consist of two echoes: one from the front surface and 
the other j'rom the back surface. The amplitude of the first one would 
give an approximate value of the dielectric constant at the surface. 

The time delay between the two echoes and the size of the nucleus would 
give the average dielectric constant throughout the nucleus, and the 
relative amplitude of the two echoes would give a gross estimate of 
the attenuation constant, 'ihese measurements, taken together, would 
give some information about the percentage of impurities and dust in 
the nucleus. 
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Figure U-8 


Figure 't 



Loss Tangent of Ice with Inpurities. The loss tangent 
of ice is plotted as a function of inpurity concentration. 
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. Fchoes from a Lossy Nucleus Signal-to-noise ratio 
of the echo from the back side of a dirty ice nucleus 
as a function of the total loss tangent. It was assumed 
that >00 echoes were added Incoherently. The distance 
from the spacecraft to the nucleus center is >0 km. 
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Tabic 4 -«. Lobs Tangent of Different Bocka* 



hock Type 

Tan ft 


Plagioclare 

3 X 10"^ 


Augitr 

3 X 10"‘ 


hypt rsthene 

o 

1 


Hasa 1 t 

o 

1 


Grunitf 

3 X 10"^ 


Dunite 

3 X 10”^ 


Liironiti 

1 

O 


Gypsum 

S X 10"** 


Orthoclase 

b X 10"' 


Olivine 

io"" 


1-luscovite 

3 X 10"** 


Palagonitc 

15 X 10"‘^ 


Pyrophylllte 

o 

1 


Chondrit es^ 

2 X 10"’ to 2 X 10"^ 

^hro.11 t. 
hew Vork, 

1. Frkhamenko, Uffitrlfial 
19b7. 

PDJD'TtUa.Of h‘ 'ks. Plenum Prt sr. 

'’Froi .1 The 

-liQCDi tdlted by Z. Kopal 

and 1 . K. F'lkhallov, Academic Pres: 


London, 19b2. 


ll the nucleus har a rocky core, the second echo would appear inuoh 
closer to the first one than in the care cl a hcno<reneour iceball. The 
tem^iorai separation between the two echoes would rive the thickness of the 
let mantle, and therefore the size of the core. The sane applies for the 
case of a roexy crust with any icy core. However, the relative and 
aosolute strength of the successive echoes will be different than in the 
previous case. This would allow us to identify the actual nature of the 
core and the crust . 

11 the nucleus has a large number of layers, then the sounding 
echoes will be similar to what is shown in Figure 2-2. 
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Because of the spherical iteonetry of the nucleus, it Is possible 
that sone Mive rays ^et refracted and returned back strairtht to the 
antenna (tiiture 4>1b). The anisles at which this effect occurs can 
be easily determined. A wavelet emitted at angle a rotates through 
angle when it refracts, and through angle n • when it reflects. 

Calling N the number of path segments within tne nucleus (e.g., the 
number of reflections is N • 1) and h the number of times the ray makes 
a full rotation, then the backward refraction condition is 


2(0, - 0^) ♦ (N - 1)(n - ?02> t ♦ 2H i (4-9) 


Also, since 

d 

sin 0 , > - sin (I 


Snell 's law gives 

d 

- sin u « sin t »2 


Combining the above equations, we obtain the condition for "back to 
the antenna" refraction: 


( d V /d \ (N-2 - 2 M)w 

- sin a 1 - N arc sin I ■■ sin a I - a ♦ . .. >0 (4-10) 

a f \av« / 2 


The zeros of this equation, which can be found numerically, 
give those values of a that will return a signal straight back toward 
the antenna. 


The lengths of the two segments of the signal path are given 
by (see Figure 4- 1b) 


= 2a cos O 2 

^out = ~ *'l 


which allows computation of the geometric attenuation. 


(4-11) 

(4-12) 
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Calculating the change In attenuation due to changes iii tiiC radii"* 
of curvature of the pulse is considerably more difficult in tnis spheric*)] 
geometry than in the previous slab model. Keeping for the present the 
approxination that the geometric attenuation Is given by the inverse 
rquare of the total path length, the returned signal as a function of time 
may be computed easily. The transmission T(0) and reflection R(0) art 
independent of the direction of the ray; that is, a ray refracts through 
the same angles and is attenuated by the same amount whether it is 
entering or leaving the nucleus. Therefore, for a given mode (i.e., a 
certain value of N), the total attenuation is 


T^rN^I 

- e 

4n(2/out ♦ N»in)^ 




where accounts for the absorption loss throughout the nucleus. R 

and T are the reflection and transmission coefficients at the surface 
(Born and Wolf, 1965). The time at which the signal arrived at the 
antenna is given by 


There are obviously several competing factors in the first expression; 
when a increases R and N increase, while T decreases. It is net obvious 
which modes give the strongest return. Computational results (below) 
show, however, that in general the returned power always decreases as R 
increases, except when the signal is very strongly attenuated. 

The computation involved in determining the refracted signal is not 
difficult. The basic approach is to iterate a, M, and N in Equation (M- 
10) and find the zeros; the values of a so .btained correspond to the 
refracted modes. Then using Equations (<*-11) through (4-14), it Id straight- 
forward to find «2, fint fout • delay t, and the total 'attenuation. 

The signal-to-noise ratio for lour main refracted modes is shown in Figure 
4-10. Curve a corresponds t*" the mode where the wave ray is along the 
line from the sensor to the nucleus center. Curve b corresponds to the 
refracted mode shown in Figure 4-1b. Curves c and d correspond to the case 
where N s 5 and 4 respectively, and M = 1. It is clear that the two first 
modes are the main contributors to the returned echo. 


D. SURFACE bACKLCATTER 

The surface of the nucleus is most probably not a perfectly smooth 
sphere. The surface roughness will backscatter some of the incident energy; 
however, because the radar wavelength is at least several meters, the surface 
backscatter contribution is expected to be very small, unless the surface is 
extremely rough. In any case, such backscattered energy will be received 
immediately after the specular echo from the front surface. Thus, such echoe 
will only affect returns from the regions immediately below the surface. 
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Figure 4-10. tchoes from the Refracted Modes. Signal-to-nolse 
ratio of four main refracted modes as a function 
of the distance between the sensor and the nucleus 
center. The nucleus is assumed to have a radius 
of 2.5 km and the loss in the ice is assumed to be 
2 dB/km. Also marked on the curves is the time of 
arrival of the echoes expressed in m** 
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SYSTEM PARAMEIiRS AND SENSOR CONFIGURATION 


Electromagnetic sounders with capabilities for deep soundini^ 

(many k'.lometers) of terrestrial ice sheets from airborne platforms 
have been built by a number of researchers. In the case of a cometarv 
nucleus sounder, the main constraints are the oower available and the 
wei(;ht limitation. 

A summary of the parameters of a point desien sounder are <<iven 
in Table M-1. These parameters correspond to a very realistic system 
compatible with spacecraft requirements and hl^th quality information 
return. The selection of these parameters was a result of the followlnr 
trade-offs: 

(1) As discussed in Section III, the optimum frequency band 
for sounding the nucleus is between 30 and 100 MHz. A 
point design frequency of 40 MHz was selected. Much higher 
frequencies were avoided because of the fact that if rocky 
debris of a fraction of a meter size are imbedded in the 
icy nucleus, they would appreciably scatter waves at hleher 
frequency, and therefore limit the soundinr depth. 

(<?) A high peak power is desirable because It would allow sorndin*' 
while the spacecraft is still at a large distance from 
the nucleus (see below). A peak power of 500 W was selected. 
Such a power can he achieved with present-day solid-state 
transmitters. 

(j) The bandwidth was selected to achieve a range resolution 

capability of better than 45 m inside the nucleus, better 
resolution can be easily achieved. 

Based on these parameters, the onorating range of the sounder 
(i.e., maximum distance between the nucleus and the spacecraft) can 
then be derived. It is reasonable to impose a requirement that the 
echo from the front surface must be 60 dll above the noise (i.e., SNR 
= bO db). Based on this requirement, the radar equation for specular 
return gives 

GA /rPp 

d s — ; — as 200 km 

dm kTB 

where r is the Fresnel reflectivity of the surface and k is the Boltzmann 
constant . 

The above point design system is capable of sounding an ice nucleus 
from any distance less than 140 km away. A higher SNR can also be 
achieved by pulse compression or incoherent addition of successive 
pulses. Evidently, for a nucleus with dirty ice, closer distance is 
required . 
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A blocK (liaRram of the point design sounder is shown in Kigure 
5-1. The sounder is essentially a pulsed radar system operating in 
the lower portion of the high frequency radio spectrum. The system 
consists of six subsystems which are as follows: 

(1) Transmitter 

(2) Receiver 

(3) Timing Control 

(4) Data Processing 

(5) Power Supply and Switching 

(6) Antenna, Matching, and TR Switching 

The transmitter transmits a pulse of known frequency at some 
pulse duration and repetition rate. After the pulse is transmitted, 
the returning echo is received and transformed to a convenient intermediate 
frequency. The signal is amplified, envelope detected, and processed 
through a data processing system and stored for transmission to Earth 
via spacecraft radio. 

The timing control circuitry insures that the transmit-recei ve 
(TR) switch is in the transmit mode when a pulse Is about to be sent, 
and in the receive mode when a returning echo appears at the antenna. 

It also controls the IF amplifier gain such that gain Increases as 
a function of time after the transmitted pulse is sent. 

The data processing system processes the video signals containing 
the echo information into a form suitable for retransmission back to 
Earth. The power supply and interface contain all of the peripheral 
components which, though necessary for proper operation of the system, 
do not belong to any one subsystem. These include items such as the 
power supplies, telemetry for systems monitoring, command telemetry 
band switching, decoupling and packaging. 

The antenna subsystem contains the antenna associated matching 
network, as well as deployment devices and the TR switch. A detailed 
description of the main subsystems follows. 


A. TRANSMITTER 

1. Master Oscillator (MO) 

The master oscillator is the single source of frequency 
generation and timing control for the entire system. It is a fixed- 
frequency, crystal-controlled, temperature-compensated type capable 
of at least 1-mW output. The unit also exhibits good aging (long-term 
stability) characteristics and frequency stability without regard to 
supply voltage or temperatures existing in a spacecraft environment. 
Commercially available units exist which would fulfill this requirement. 
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Figure 5-1. Block Diagram of the 


Point Design Sounder 
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Buffer Amp 


The oscilla^or buffer amplifier provides the master oscillator 
with isolation from the rest of the system, Insurlnc; a stable master 
frequency source, free from the effects of frequency "pullln*?" which 
could occur durlnq periods of peak power supply demands. This unit 
also provides a nominal power gain to offset the effect of driving 
the digital divider chain and local oscillator mixer. 


3. Second Buffer Amp 

The second buffer amp amplifies the output of the first buffer 
amp to a level sufficient to drive the local oscillator mixer (7 dB). 


4. Isolation Amp 

The Isolation amplifier amplifies the output of the first hut fer 
amp (or frequency divider) to a level necessary to drive the puise 
modulator. 


5. Digital Divider Chain (OPD) 

The digital divider chain Is clocked by a portion of the first 
buffer amp output. It divides the 20 MHz clock by 2 x 10*’ resulting 
In an output of 100 pulses/s. 


b. Pulse Repetition Rate Generator (PRR) 

The pulse repetition rate genei ator Is a divider (+10) which, 
when enabled by a command pulse, fires the pulse width generator. 

It also provides timing pulses to the TR switch tlmlug control and 
AGC timing control alerting those units that a pulse is about to bo 
sent or received. 


!. Pulse Width vJenerator (PWCJ 

The pulse width generator clocks the pulse modulator for a precise 
length of time allowing it to transmit the carrier frequency to the 
pulse modulator input. 


d. Pulse Modulator (PM) 

The pulse modulator is a double-balanced mixer (DhM) which is 
utilized as a low-level modulator. It accepts the TTL digital level 
source from the PWd and carrier frequency from the ISO amp. It outputs 
t)io carrier frequency to drive the first RF driver -asn only when a pulse 
is available at the mixer input. The output is anproxlmately -3 dP. 
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9. Hrst Driver Anp (DAI) 

Ine first driver .unp is sourced from th€' Ff'. output. It provides 
about 25 db ftain, giving an output of approximately 22 dbm which drives 
the second driver amp. 


10. Second Driver Amp (DA2) 

The second driver amp i.s a class ”C" ;«npllfier with 25 dF >^iln. Class 
*'C" opcratlor gives n mueii hlglu r powt-r eificiency since it Is l>la.«d ol f 
during quiescent opt rat ion. The p<jwer output is approximately 50-lw peak. 


11. hi Fowtr Amp (FA) 

The HF power amp is the final power amplltler. It is a class 
"C" Diased amplifier requiring 50-w drive and having an output of 500-W 
peak (gain x 10 db) . Solid-state design is ii.srd to eliminate t lie 
requirement for high-voltage power supplies. The output of the amplllier 
is fed directly to the TR switch to the antenna. 


b. SbCtlVbh (hX) 

1. Input (•liter (I PI) 

The input filter is located immediately after the TM switch, 
it receives the broadband Hi noise and signal from the antenna. The 
bandwidth is sufficient to pass t(ie des ired returning echo of tne transmitted 
pulse and reject the unwanted Image frequency and broadband noise. T(ie 
output directly feeds the Hi amplifier. Addition of tliis filler increan i 
Uie absolute sensitivity of tbo system but also adds directly to the 
noise figure of ♦he system. The loss is expected to be lesr i han I.*! db. 


2. Hi Ampliiier (Hil) 

The Hi amplifier i.s a solid-state broadband design with wi le 
dynamic range (80 db). The nominal gain is approximately 14 db and 
noise figure of 4 dP. 


3. First Mixer (MX1) 

Tne first mixer Is a double-balanced broadband type which pi-ovides 
superior isolation. Tlie Hi output of the Hil ar.pl I tier i.n mixed with 
the local oscillator and the difference pn duet t(ii mixer output 
is selected as the IF frequency. 



4. IF Filtera (FF) 

There ia an IF filter at the carrier frequency with an 8 MHz 
bandwidth. 


5. IK Amplifier (IFA) 

The IF amplifier providea 80 dB amplification of the IF aignal 
to the video. Proviaiona for 60 dB of gain control are included. 
Amplifier bandwidth (3 dB) ia 10 and IF frequency ia 30 MHz. A 
video detector and video amplifier are alao incorporated into the deaign. 
Video amp gain ia 10 dB giving a total of 90 dB aignal gain from the 
IF input. 


6. AGC Timing Control (AGC) 

AGC timing control ia provided to vary the IF gain aa a function 
of time ao that the video level will remain conatai.’ with echo range. 


7. Second Mixer (MX2) 

The aecond mixer ia a double-balanced type. Thia mixer adda 
the local oacillator frequency to the carrier frequency to generate 
a mixing frequency for MX1. Ita high iaolation providea 30 dB auppreaalon 
of the local oacillator aignal. 


8. Local Oacillator (LO) 

The local oacillator is a stable crystal oscillator type with 
sufficient output to drive MX2 (approximately equals 7 dBm). The frequency 
is tht same as the IF frequency. 


9. LO Bandpass Filter (LOF) 

The local oscillator bandpass filter will have narrow selectivity 
and high shape factor. This is to insure that no unwanted mixer output 
product passes. Unwanted products are expected to down 30 dB through 
this device. 


10. Local Oscillator Amp (LOA) 

The local oscillator amp amplifies the filtered output of MX2 
to provide sufficient power to drive MXI. 
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The transmlt-receive switch selects either the transmitter output or 
receiver Input to the antenna. Isolation should be 60 dB to provide suffi- 
cient protection to the receiver input and allow RFA to recover quickly 
enough to receive the returning echo. Switching time should be to 100 ns. 


2. Switch Timing Control (STC) 

The switch timing control is used to control the TR. It must 
activate the switch fast enough to Insure that no part of the transmitted 
pulse enters the receiver directly from the transmitter. 

3. Antenna Matching Network (AMN) 

The antenna matching network transfers the maximum amount of 
power to the antenna through a tuning network and Baiun transformers. 
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CONCLUSION 


The potential significance of science returns from deep sounding 
of a cometary nucleus is very high. Rlectromagnetic echoes returned 
from the comet can be interpreted to deduce the internal structure 
of the nucleus. The cometary plasma places a lower bound of 10 MHz 
on the sounding frequency; ice absorption requires the frequency to 
be below 100 Mhz for deep penetration. 

Hadiatlon scatters back to the antenna by a variety of mechanisms; 
the most important are reflection from layered dielectric interfaces 
and Inhomogeneities, and radiation refracted through the nucleus. The 
strength of these signals and their time delays yield a profile of 
the dielectric changes throughout the nucleus. 

A point design sounder which can operate from a 1U0>km distance 
at 40 Mhz is discussed in some detail and found to be within the current 
state-of-the-art of space technology. 

In conclusion, a sounder can be developed for the payload of 
a rendezvous comet mission, and its scientific payoff is appreciable. 

It is basically the only sensor which can accurately determine the 
complete internal structure of the nucleus. 
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